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I.  INTRODUCTION 


The  Spectral  Theory  of  Diffraction  (STD) ,  which  was  introduced  recently 
by  Mittra,  Rahmat-Samii  and  Ko  [1],  has  been  found  to  be  a  useful  tool  for 
solving  high-frequency  diffraction  problems.  The  concepts  of  STD  together 
with  its  scope,  advantages  and  disadvantages  can  be  found  in  a  number  of 
publications  [1],  [2]  which  have  appeared  recently  in  the  literature.  This 
theory  is  founded  on  the  concept  that  the  scattered  field  is  related  to 
the  Fourier  transform  of  the  induced  current  on  the  scatterer  and  that 
the  ray  description  of  the  scattered  field  is  derivable  from  the  spectral 
representation  via  an  asymptotic  evaluation  of  the  Integral  form  of  this 
representation. 

The  canonical  geometries  of  the  semi- inf inite  half -plane  and  the  wedge 
play  important  roles  in  the  development  of  the  various  theories  for  high- 
frequency  scattering.  The  Sommerfeld  half -plane  solution  has  been  extensively 
studied  [2]  in  the  context  of  STD,  and  the  STD  formalism  for  the  half-plane 
diffraction  has  been  applied  to  planar  and  curved  strips.  In  this  paper, 
we  extend  the  STD  concepts  to  the  geometry  of  the  wedge.  We  show  how  the 
contour  integral  a  la  Sommerfeld  for  the  wedge  problem  can  be  manipulated 
to  yield  a  transform  representation  for  the  scattered  field  in  terms  of 
the  currents  flowing  on  the  two  faces  of  the  wedge.  We  also  show  that  it 
is  possible  to  derive  a  uniform  asymptotic  solution  for  the  total  field, 
one  which  is  valid  for  all  observation  angles.  The  uniform  solution  so 
derived  is  found  to  be  identical  to  that  conjectured  by  Deschamps  and 
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The  starting  point  of  our  analysis  is  the  Sommerfeld  solution  to  the 

problem  of  plane-vave  diffraction  by  a  wedge.  Let  a  perfectly  conducting 

wedge,  with  the  faces  defined  by  the  angles  $  ■  0  and  $  ■  $,  be  illuminated 
* 

by  a  plane  wave 

ui  -  e-ikPcos<+-*')  (1) 

incident  from  the  <j>*  direction.  (See  Fig.  1.) 


Figure  1.  Geometry  of  the  wedge-diffraction  problem. 

The  angles  involved  have  the  following  ranges: 

3 

<  2ir 5  0  £  $  £  0  <  <j>’  <  ir. 


The  total  field  u  is  subject  to  the  boundary  condition 

u  -  0  on  the  wedge  surface,  for  u  ■  E  (E-wave) ,  or 

z 

*  0  on  the  wedge  surface,  for  u  ■  H  (H-wave). 
dn  z 


*  —  Itot 

e  time  convention  is  implicit  throughout. 
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This  field,  as  derived  by  Sommerfeld,  can  be  represented  in  the  form  [3] 


where 


u(p,<j>)  -  w(p ,<J> ; * )  +  w(p,<|> 


w(p,$;<j>) 


_2_  | •  -ikpcos(a-4>)  da _ 

2it  aJ  1_ei£i(a-ii») 


(3) 

(4) 


is  the  Sommerfeld  integral  defined  on  the  contour  A  in  the  complex  a 
plane  (Fig.  2),  and  -  j.  The  upper  (lower)  sign  in  (3)  applies  to  the 
case  of  the  E-wave  (H-wave) . 


Figure  2.  Contour  A  ■  in  the  complex  a 

plane.  The  shadowed  areas  indicate  the 
convergence  regions  for  the  Sommerfeld 
integral  w(p,$;$). 

III.  ALTERNATE  FORM  OF  THE  SOMMERFELD  INTEGRAL  w 
In  this  section  we  derive  an  alternate  form  of  the  Sommerfeld  integral 


w,  one  which  would  allow  us  to  express  the  scattered  field  in  terms  of  the 
transform  of  the  induced  current  on  the  two  faces  of  the  wedge. 
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We  begin  with  the  integrand  for  w  given  in  (4)  and  note  that  the  poles 
are  located  on  the  real  axis  in  the  a  plane  at  the  points 

otp  ■  ip  +  2p$,  p  •  0,  ±  1,  ±  2,  . . . 

Since  <t  >  n,  only  one  of  these  poles  can  appear  in  the  region  (— -rr  + 
n  +  <p) .  This  implies  that  only  the  incident  and  one  reflected  wave  can  reach 
the  given  observation  point  (p,<J>).  The  next  step  is  to  replace  the  contour 
A  by  a  sum  of  the  contours  K  and  L  (■  +  I^) ,  as  shown  in  Fig.  3. 


Rea 


Figure  3.  Contours  K  and  L  (=  +  I^)  in  the  a  plane. 


Because  the  contours  and  are  separated  by  2ir,  and  traverse  in  opposite 
directions,  we  can  write 


Jf(a)da  *  -/f  (a+2ir)da. 


Consequently,  the  integral  along  L  can  be  replaced  by  one  along  either 
or  L2‘  Next,  let  us  turn  to  the  integral  along  the  contour  K.  The  contri¬ 
bution  to  w  from  the  integral  along  K  may  be  found  by  evaluating  the  residue 
at  the  pole  a^.  Incorporating  these  manipulations  in  the  expression  for  w 
given  in  (4) ,  we  get 


(5) 


w(p,<j>;<iO  -  En(ir-[ap-ij>])e  ikpcos(ap  +  l(p , 


where 


r _ i_  t  -ikpcos(ct-<i>)  sinflir  da 

I  2u  ^  poaflfn— i 


cosO  (a-ijrHr)  -cosflir 


.  or 


i(p -  I 


1  j  ^-ikpcos(a-<|))  0  ainflir  da 


^  2lT  L 


cosQ(a-^-ir)cosftir  * 


(6a) 


(6b) 


and  n(x)  Is  the  Heaviside  step  function.  The  sum  in  (5)  is  extended  over 

all  p’s  for  which  the  inequality  it— £oi^— d> ] >0  holds  at  least  for  some  value 

of  <|>  lying  in  the  range  0  <  <j>  <  $. 

Simple  asymptotic  considerations  reveal  that  the  sum  over  p  in  the 

r.h.s.  in  (5)  is  of  the  order  of  (kp)°  as  kp-*»,  whereas  the  integral  I  is 

-1/2 

of  the  order  of  (kp)  .  Hence  the  representation  (5)  can  be  useful  for 
high-frequency  asymptotic  analysis.  Also,  it  will  be  seen  from  the  following 
that  the  integral  I  can  be  easily  converted  into  Fourier  Integral  form. 

Both  of  these  features  will  be  taken  advantage  of  in  subsequent  sections. 


IV.  RAY-OPTICAL  REPRESENTATION  OF  U 

When  <f>*  and  <j>  lie  in  the  range  0  <  <j>’  <  ir,  0  <_■$<_  #,  only  one  pole 
a  ■  aj  of  w(p,$;$’)  and  two  poles  a  ■  <Xq,  and  a  »  of  w(p,<J>  ;-<!>')  can 
satisfy  the  inequality  -n+4><a<ir+$ .  These  are 

<Xq  ■  due  to  the  incident  wave, 

Oq  ■  -<p' ,  due  to  the  wave  reflected  from  the  wall  <j>  *  0,  and 
of  ■  2$  -  <p' ,  due  to  the  wave  reflected  from  the  wall  d>  ■  *. 
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Thus,  we  find  from  (5) 


w(p, $;-$') 


W<p , 4> ; 4> * )  *  n(n+4>,-<)>)e"ikpCOS^  ^  »  and  (7) 

-  n[«-(*'+*)]e'1‘“!cos<*''h>)  +  ^2J-(*'++)-.]rll®cos<M-*'-'») 


Substituting  (7)  and  (8)  into  (3) ,  we  obtain  the  representation  for  u 

which  is  amenable  to  ray-optical  interpretation.  Indeed,  the  first  term 

in  (7)  and  first  two  terms  in  (8)  are  of  the  order  of  (kp)°  as  kp-n»,  and 

represent  the  GO  approximation  to  the  total  field.  The  remainders  in  (7) 

and  (8),  represented  by  I(4>’)  and  I d> 1 )  ,  respectively,  are  of  the  order 
-1/2 

of  (kp)  and  represent  the  field  diffracted  at  the  edge  of  the  wedge  (for 
details  see  [4],  [5]). 


V.  CONVENTIONAL  REPRESENTATION  OF  U 

The  representation  for  the  total  field  u,  given  in  (3),  (7)  and  (8), 
contains  terms  which  are  discontinuous  functions  of  the  observation  angle 
(j>  though  the  total  field  is  continuous.  It  is  useful  to  seek  an  alternate 
representation  in  the  context  of  STD,  which  consists  of  continuous  functions 
only.  This  is  accomplished  by  deforming  the  contour  L  such  that  the  poles 
ctg,  a.Q  and  do  not  cross  the  new  contour  as  and  <j»*  are  varied.  We 
choose  this  deformed  contour  to  be  for  $  <  it  and  L*  for  <J>  >  ir.  Both  of 


these  contours  are  shown  in  Fig.  4. 
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VJ 


0  <p 


1 

QS 

j 

Figure  4.  Contours  L  (<j><ir)  and  L"  (<j»>ir)  in  the  a  plane. 


With  such  a  choice  for  the  contours,  the  pole  is  always  located  inside 
the  region  bounded  either  by  L*  and  L^,  or  by  L*  and  L^,  and  the  poles 
and  a*  always  appear  outside  of  the  regions. 

The  total  field  u  can  now  be  expressed  as 
i  ,  s 

i  s 

where  u  and  u  are,  respectively,  the  incident  and  scattered  fields.  The 
scattered  field  is  defined  by  the  formula 


s  iA  f  -ikpcos  (a-<j>)17,  .  . 

u  "  J  e  W(cO  da» 


where 


\  -  sgn  y,  P. 


A  -  1 


L2  A  -  1 


and  W(a)  is  defined  as 

W(u)  ■  OsinfliK [cosfi(a-t|H-ATr)-cosn7r3  ^  +  [cosft(a+i|rK\Tr)-cosfi7r]  (1 


the  upper  (lower)  sign  corresponds  to  the  case  of  the  E-wave  (H-wave) . 
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Our  choice  of  and  for  the  contour  ^  is  arbitrary  in  that  we 
are  free  to  use  either  or  for  P_  and  or  for  P+,  with  appro¬ 
priate  choice  for  the  corresponding  integrals. 

In  the  literature,  the  expression  given  in  (9)  is  referred  to  as  the 
conventional  representation.  (See  reference  [1].)  We  proceed  next  to 
find  the  Fourier  transform  representation  of  the  scattered  field.  This 
is  described  in  the  following  section. 

VI.  FOURIER  TRANSFORM  REPRESENTATION  OF  THE  SCATTERED  FIELD  u8 


Let  us  introduce  the  following  change  of  variables: 

/~2  2 

w»kcosa;  v*/k  -w  *-Xksina;  x**pcos<|>;  y-psimfc; 
k  “-kcos<j>  ’ ;  k  --ksin<j> ' . 

^  y 

For  v,  the  branch  of  the  square  root  is  chosen  such  that  Im  v>0. 
Since  the  transformation  w  a  is  a  multiple-valued  function,  the  two 
contours  and  in  the  a-plane  transform  into  a  single  contour  C  in 
the  w-plane.  The  resulting  contour  in  the  w-plane  is  shown  in  Fig.  5. 


Rew 


Figure  5.  Contour  C  in  the  complex  w-plane  (image  of 
L*  and  in  a-plane) . 
w.  ■  -k  ,  w.  ■  -(cos2$k  +  sin2$k  ) 


Thus  the  scattered  field  uS,  as  defined  by  (10)  and  (11),  becomes 


If  we  define  the  direct  and  inverse  Fourier  transforms  for  the  functions 


u(x)  and  U(w)  as 


U(w) 


/  e*WX  u(x)  dx; 


u(x) 


,  ®+iA 

i  / 


-OO+iA 


e  *wx  U(w)  dw,  A>0, 


(15) 


then  a  comparison  of  (13)  and  (15)  reveals  that  uS  in  (13)  has  the  form 

eiv I y I 

of  an  inverse  Fourier  transform  for  the  function  — — —  2V.  We  observe 

.iv  M  21v 

that  — 2iv~  is  the  Fourier  transform  of  the  free-space  Green's  function 
^  (kR)  and  that  the  scattered  field  has  been  expressed  in  the  form  of 
a  convolution  Integral  involving  the  Green's  function  g  and  the  surface 
current  j .  These  integrals  take  the  form 


m 

u  (r)  »  -lain  /  g(£-r')j (r')dr'  for  E-wave,  and  (16a 

s 

u3(r)  ■  ■—  /  g(r-r^')j  (r')dr'  for  H-wave.  (16b 

s 

Returning  now  to  (13)  we  deduce  that  the  function  V(w)  is  proportional  to 
the  Fourier  transform  of  the  current  flowing  on  the  surface  of  the  wedge. 
This  conjecture  will  now  be  verified. 


VII.  ALTERNATE  FORMS  FOR  THE  SCATTERED  FIELD 

In  this  section  we  derive  alternate  forms  for  the  scattered  field  which 
can  be  interpreted  in  the  context  of  the  spectral  domain.  We  consider  the 
E-wave  and  H-wave  cases  separately. 


7.1.  E-Wave 


From  (16)  It  follows  that  the  scattered  field  u  can  be  represented  as 
a  summation  of  two  terms,  viz., 

u  -  uQ  +  (17 


where 


Ug(r)  -  -low  /  dx  j  Hq1^ (k|r-fcx| )jQ(x) , 


00 

u®(r)  =  -lay  /  cb^  ~  Hg1^  (k|r-*1x1j)j1(x1), 


The  tiro  terms  In  the  r.h.s.  of  (17)  may  be  associated  with  the  two  currents 


JnW  “  “  7Z7 


1  3u< 


ia>li  3y 


1  3»(x1,y1-0) 

"  la, 

which  flow  on  the  two  surfaces  of  the  wedge  defined  by  <p  ■  0  and  <(>»$, 
respectively  (Fig.  6). 


GOG  ]0W 


Figure  6.  The  currents  jg  (x)  and  on  the 

wedge  surfaces  defined  by  <p  ■  0  and  <(>  *  4>, 


respectively.  Both  currents  flow  along  the  z-axis. 


Replacing  the  Hankel  function  In  (18)  by  Its  Integral  representation  [4] 


i  ---hi 


^  ^ivKx-x'J+ivly-y' 


2  2 

dw;  v  *  k  -w  ;  Im  v>0; 


^x-x')2  +  (y-y*)2 


and  changing  the  order  of  Integration  leads  to  the  following  expressions  for 


s  s 

Uq  and  u^: 


J*\y 


v->-- 


(20a) 


s  /  \  ium  r  ,  -lw,x,  e^Vl^l  „  ,  ^ 

Ul(^  "  -  17  l  dwl  e  11  "  '  2iv, J1(W1>  *  Where 


(20b) 


Jo  *  /  dxe  ^o  ^  *  and 


(21a) 


Jl(wl)  -  /  dx-je  wlXlj1(x1) 


(21b) 


3  3 

We  observe  from  (20)  that  the  component  fields  u^  and  comprising  the 

3 

scattered  field  u  have  the  form  of  inverse  Fourier  transforms  of  functions 
related  to  the  Fourier  transforms  of  currents  jg  and  j  ^  flowing  along  the 
faces  of  the  wedge.  This  is  one  of  the  results  we  had  been  seeking. 

Next,  we  show  that  the  scattered  field  u  can  be  written  as  a  Fourier 
transform  of  a  single  quantity  which  is  related  to  a  suitable  combination 
of  the  two  currents.  We  begin  by  introducing  an  additional  change  of 
variables,  defined  by 


w,  -  k  cos(a-$),  -X,v,  -  k  sin(a-$),  where  X,  -  sgn  y, . 
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Such  a  change  of  variables  Is  equivalent  to  the  enforcement  of  the  condition 
-w^  +  vjyj  -  -wx  +  v | y |  . 

This  may  be  verified  by  using  the  relationship 

x^  -  xcos$  +  ysin$,  y^  -  -xsin$  +  ycos4>. 

When  these  changes  of  variables  given  In  (22)  are  inserted  Into  (20a)  we  get 


4-  I  <*«  ,-:U«c08ta-'*)  Jg(kcosa) , 


(23) 


where  P.  is  defined  in  (10a) ,  and 


X-u>U 


r 


U1  * - h~  /  do  e  ^,^pco8^ot  J^[kcos(a-(^) ] ,  where  L  -  / 


L3  Xj-1,  X-l, 

14  xx— 1,  X-l,  . 

15  xx— 1,  X— 1, 


The  contours  L3>  and  L5  are  shown  In  Fig.  7, 
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The  contour  can  be  deformed  into  the  contour  (see  Fig.  4) ,  the 

contour  L,  into  the  contour  which  coincides  with  L,  but  is  orientated 
4  l 

along  the  opposite  direction,  and  the  contour  into  the  contour 
With  these  changes,  u®  can  be  rewritten  as 

u?  -  -  ^  /  da  .-itocos(o-+>J,[kc0a<a-*>). 

1  ZTT  p  1 

Using  (14)  and  (23)  in  (24),  we  finally  obtain 
us  ■  -  /  do  e  ikpcos(a  lip.  f  (Iccosct)  +  J^[kcos(a-<&)]  ] }. 


Comparison  of  this  result  and  (10)  gives 


W(a)  *  -  — p-  (J  (kcosa)  +  J^fkcosfa-*) ] } •  (26) 

This  formula  is  useful  for  providing  a  physical  interpretation  of  the 
function  W(a)  in  that  the  scattered  field  uS  is  a  Fourier  transform  of  W(o) 
which  is  related  to  the  transforms  of  Jq  and  as  shown  in  (26) .  Note 
that  equating  the  integrands  of  (10)  and  (23)  is  permissible  because  the 
two  integral  representations  for  uS  are  valid  for  any  contour  resulting 
from  the  deformation  of  P^  in  Fig.  5. 

Next  we  proceed  to  separate  the  physical  optics  (PO)  and  the  fringe 
currents  in  the  transform  domain.  We  begin  by  introducing  the  function 
Wf(o)  which  is  to  be  interpreted  later  as  the  transform  of  the  fringe  current: 

(a)  »  W(a)  -  W^(a),  (27) 

where  WP0(a)  is  obtained  from  (26)  by  replacing  Jq  and  with  and 
PO 

,  respectively,  the  latter  being  the  transforms  of  the  physical  optics 
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.5 


one  can  show  that  the  singularities  in  (27)  cancel  each  other  out  and, 

consequently,  W^a),  the  transform  of  the  fringe  current,  is  a  regular 

function.  If  we  replace  W(a)  in  (10)  with  (a)  and  apply  the  method  of 

stationary  phase  to  the  resulting  integral,  we  obtain  a  field  of  the  order 
-1/2 

of  (kp)  which  is  regular  away  from  the  wedge,  and  includes  the  shadow 
boundaries  of  the  incident  and  both  reflected  waves.  This  field  can  be 
attributed  to  the  fringe  current,  which  is  the  excess  over  the  PO  current 
on  the  surface  of  the  wedge. 


7.2  H-wave. 

The  scattered  field  for  the  H-polarization  case  can  be  again  expressed 

s  s 

by  (17),  but  the  function  u^  and  u^  now  becomes 

00 

and 

ui(£)  ■  ■  3?7  J  dli  I  Bo1)<kl£-£iltil) 

where 


jQ(x)  «  u(x,y  -  0)  and  j^x^  -  u(x1,y1  -  0). 

The  current  jQ  (x)  flows  along  the  x-axis  and  the  current  j^(x^)  along  the 
negative  x^  -  axis.  (Fig.  8) 

By  applying  a  procedure  similar  to  the  one  for  the  E-wave  case,  we  find 


u8  ■  -  — >  /  da  e  iltPcos(a  {JL  [8ina  Jn  (kcosa)  -  sin(a-$)  J.  [kcos(a-<P)]  ] }. 


(29) 
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Figure  8.  The  currents  jQ(x)  and  j  ^  (ac^)  flowing  on 
the  surfaces  $  ■  0  and  $  ■  #. 


Then  from  (10)  and  (29)  we  obtain 

W(a)  -  {sinaJ^j  (kcosa)  -  sin(a-t) J^kcosfa-#) ] } . 

The  function  W(a)  can  again  be  interpreted  as  a  combination  of  terms  related 
to  the  Fourier  transform  of  the  currents  on  the  faces  of  the  wedge.  All 
further  observations  regarding  W* (a)  continue  to  remain  valid,  except  that 
the  expressions  for  and  J?°  are  now  given  by 


J()0 (kcosa)  "  cosa-cos*’  ’  *®<cosa-co8*')>0  ,  and 
J^°[kcos(a-$)]  *  Y  Cos(a-«) -cos ($'-*)'  ’  I®  [  cos  (a-*) -cos  (<►’-*)  ]>0. 


The  main  result  of  this  section  is  that  we  represented  the  scattered  field 

g 

u  in  terms  of  Fourier  transforms  of  currents  flowing  along  the  faces  of  the 
wedge.  Even  though  the  transforms  are  discontinuous  functions  in  the 
spectral  domain  (compare  (26)  and  (11)),  the  integral  representation  for  the 
scattered  field  (25)  remains  a  uniformly  continuous  function  of  the  observa- 
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tion  angle  <j>  In  Its  range  of  variation.  Of  course,  the  total  field  is  also 
a  continuous  function  of  this  angle. 

In  the  next  section  we  construct  an  approximate,  uniformly  continuous, 
asymptotic  representation  for  the  total  field  -  one  which  is  expressed  in 
terms  of  known  functions. 

VIII.  UNIFORM  ASYMPTOTIC  EXPRESSION  FOR  THE  TOTAL  FIELD  U 

As  a  first  step  toward  constructing  a  uniform  asymptotic  expression 
(UAE) ,  we  return  to  the  representation  given  in  (3)  for  the  total  field 
u(p,<t>).  Next,  we  manipulate  the  expressions  for  the  functions  w(p,<j>;<t>') 
and  w(p, $,-$')  to  derive  a  representation  which  remains  valid  for  all 
observation  angles.  We  begin  by  considering  the  function  w(p,4»4')<  For 
the  purpose  of  later  manipulation,  it  is  convenient  to  represent  this 
function  in  the  form 

w(p , 4» ; <|> * )  *  n(n-N>’-<(»)e  lkpcos^  +  I^(p ,4, ; 4> 1 )  +  I2 (p , d> ; (j) * )  ,  where  (30) 

-ikpcos(a-<f>) 

I  <P.W)  -  27  f  in(a-*')  d°* 

L  1-e 
n 

In  order  to  find  its  UAE  we  will  employ  a  technique,  described  in  [4], 
which  is  suitable  for  an  asymptotic  analysis  of  integrals  which  contain  a 
simple  pole  in  the  neighborhood  of  the  saddle  point.  According  to  this 
technique,  the  integral 

I(X)  -  /  f(a)eX8(ot)da  , 

SDP 


with  a  simple  pole  at  o  ■  and  nearby  saddle  point  at  a  *  ag,  has  the 
asymptotic  expression: 


2 

I(A)^eXg(°s'{±i2a/irVXb  Q(±ibvT)  +  «/£*((>)}  »  M^O, 


as  A  -*■  »  ,  where 


a  -  ltm[(ct-an)f(a)];  b  -  /g(a  )-g(a  )  ;  T(0)  -  hf  (a  )  +  -  ; 

a  -*■  aQ  0  3  a  8  b 


g"(ag)  * 


Q(y)  -  / 


-x 

e  dx, 


and  the  sign  of  H  Is  defined  such  that  arg  h  •  (arg  da)  ,  where  da  Is  an 

a 

s 

element  along  the  steepest  descent  path  SDP,  while  arg  b  Is  defined  such 


a_-a 

,  0  s 

b  >  — - —  ,  as  a.  +  a  . 
h  us 


In  our  case  we  have 


g(a)  -  -icos(a-<|»);  f(a)  -  ^ ;  aQ  -  4> f  >  and  on  1^: 

1-e 


a  *  -x+4;  a  - 


-  -Hr-;  h  -  e  *4/2  :  b  -  /2  e*4  cos  ; 


Let  us  assume  that  the  pole  a  *  <j>  lies  near  the  contour  L^.  We  observe  that 
i—  -iT  ®  j  2 

Q(y)  -  /rTF(e  4y) ,  where  F(z)  *  — -  /  e  y  dp  is  the  Fresnel  integral,  and 


that  F(z)  “  /iri  -  F(-z).  Consequently, 


±  b), 

and  the  sum  of  the  first  two  terms  in  (30)  have  the  following  UAE: 
n(ir+b'-««-lb,cos(*'''»)  +  I1(P,bi*0^ltetb'llWtl+COS<*''W,F(-^'<:<,81Yt'  + 


2/2irkp  cos 


-i I 

fle  4  _ 

/2irkp  1-e 


— ± -  \ 

in^-ir-*')  s  * 


Since  the  pole  a  ■  aQ  is  away  from  the  contour  L2,  the  contribution  to 

w(p ,  4» ;  1 )  from  (p ,  4> ;  d>  * )  can  be  calculated  by  applying  the  standard  saddle- 

point  method,  with  the  result 


e1^-1!  0 


>2irk p  1-e 


ifl  ($-Hr-<£ ' ) 


Upon  adding  (31)  and  (32)  together,  we  obtain  the  needed  UAE  for  the  function 
w(p ,  <|> ;  4>  * ) : 


b.  e1’“{e-lke[1+c<,,(*'-*)1F(-/2i5'cos^)  + 


il 

e  4  1 

2/2wkp  cos  ^ 


i- 

e  4 _ nsinflir _ 

/2trkp"  cosfi ($-$')  -  cosflir 


If  the  pole  ctg  lies  near  the  contour  L 2»  then  the  above  reasoning  can  be 
repeated  with  the  contours  and  L 2  interchanged.  It  can  be  shown  that 
the  final  result  is  the  same. 

To  find  the  UAE  for  the  function  w(p  ' ) ,  let  us  represent  it  in 
the  form: 


I-lCp  ,4>  ;-<!>  * )  +  i2(p. 


(33) 


fire 


where  I _  (p ,  )  **  ~Z~  / 


-ikpco8(a-(j») 


2nLJ  1_ein(a+(j),) 
n 


da 


In  this  case  the  pole  a  -  -<t>'  appears  near  the  contour  and  another 
pole  a  *  2$-<f> '  lies  near  the  contour  .  Therefore,  now  we  employ  the  same 
technique  to  the  first  and  the  third  and  to  the  second  and  the  fourth  terms 
in  (33),  respectively,  as  was  used  to  find  (32).  The  result  is 

w(p,*i-*’)^lkpte-llt(,(1+cos<*+*')1F(-/2E-cos  - ijrrr  + 

2  2/2^  cos  ^ 


-i JL 

Ue  4 


+  e 


-ikp[l+cos(2*-4«-*)]p(.^-coa  +  (34) 


e!_  1  fle^t  1  , 

2/2^  COB*t±^±  l-em"++-*'} 

The  sum  of  (33)  and  (34)  yields  the  UAE  for  the  total  field  u.  It  can  be 


rewritten  as  follows: 
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=  -pcos(ij>-<f) ' ) ;  Srl  =  -pcos (<$+({>  ' ) ;  sr2  =  -pcos(2^-<J>,-<{>) . 


In  the  context  of  ray  optics,  these  quantities  are  the  eiconals  of  the 

diffracted,  incident,  and  the  two  reflected  waves,  respectively.  We  also 

i  r  r 

introduce  the  functions  e  ,  and  which  are  equal  to  -1  in  the  lit 
region  and  +1  in  the  shadow  region.  Next  we  use  the  notation 


f\, 

F(x) 


i(x2+|-) 

e _ 4 _ 

2/tT  x 


which  is  the  first  term  of  the  asymptotic  expansion  of  F(x),  according  to 


F (x)  =  n(-x)  +  F(x)  +0  (x  "*)  . 

Then  the  UAE  for  the  total  field  u  takes  the  form: 

u  n,  ui[F(eiv4(s-si)  )  -  Jf(eiv/k(s-si)  )]+{!-►  r^}  +  {i  -*■  r +  ud,  (35) 
where 

ikp+i^-  ,  _ 

u  ,  * - -  ■  — —  2sinfiir{  [cos5H<j>-b')-cosnr]  +  [cosfi(4>+i}> ' ) -cos^tt ]  }.  (36) 

/27rkp 

We  note  that  the  expression  for  the  total  field  u  contains  u^,  given  in 
(36)  which  is  the  diffracted  field  &  la  Keller.  This  diffracted  field  is 
non-uniform  in  nature  as  it  is  singular  both  at  the  shadow  boundary  and  in 
the  neighborhood  of  p  *  0.  However,  the  total  field  given  by  (35)  is  uniform 
at  the  shadow  boundaries  because  the  singular ies  in  ^'s  appearing  in  the 
first  three  terms  exactly  cancel  out  the  singularity  in  u^.  However,  such  a 
cancellation  does  not  occur  at  p  =  u.  Hence  the  expression  for  u  in  (35) 
is  still  non-uniform  at  the  tip  of  the  wedge;  nonetheless, the  uniform  nature 
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is  maintained  at  asymptotic  distances  away  from  the  tip,  i.e.,  for  kp  large. 

Finally,  we  point  out  that  the  expression  in  (35)  is  identical  to  the 
one  based  on  the  Deschamps  and  Lee  ansatz  [6]  which  is  a  generalization  of 
the  uniform  formula  for  the  half -plane  derived  by  Ahluwalia,  Boersma  and 
Lewis  [ 7 ] . 
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